Neutron Oscillations to Parallel World: Earlier End to the Cosmic Ray Spectrum? 



O 

in 

(N 

^ : 

Oh- 

I ■ 

o ■ 

(N ■ 
> ■ 

in ; 

i> ■ 
m ; 

0^ " 

o : 



Zurab Berezhiani^' ^ and Askhat Gazizov^ 

^ Dipartimento di Fisica, Universita dell'Aquila, Via Vetoio, 67100 Coppito, L'Aquila, Italy 
^INFN, Laboratori Nazionali Gran Sasso, 67100 Assergi, L'Aquila, Italy 
^DESY Zeuthen, Platanenallee 6, D-157 Zeuthen, Germany 
(Dated: September 22, 2011) 

Present experimental data do not exclude fast oscillation of the neutron n to its degenerate 
twin from a hypothetical parallel sector, the so called mirror neutron n' . We show that this effect 
brings to a remarkable modifications of the ultrahigh-energy cosmic ray spectrum testable by the 
present Pierre Auger Observatory (PAO) and Telescope Array (TA) detector, and the future JEM- 
EUSO experiment. In particular, the baryon non- conservation during UHECR propagation at large 
cosmological distances shifts the beginning of the GZK cutoff to lower energies, while in presence 
of mirror sources it may enhance the spectrum a.t E > 100 EeV. As a consequence, a significant 
reduction of the expected diffuse cosmogenic neutrino flux is predicted. 



There may exist a hidden parallel sector that is an 
exact copy of our particle sector so that each ordinary 
particle: electron e, proton p, neutron n etc. has its mass 
degenerate twin: e' , p' , n' etc. These so called mirror 
particles are sterile to our strong and electroweak inter- 
actions SU{3) X SU{2) X C/(l) but have their own gauge 
interactions SU{3y x SU{2y x 1/(1)' with exactly the 
same couplings. Mirror baryons are viable as asymmet- 
ric dark matter provided the parallel sector has smaller 
temperature than the ordinary one T' <^ T 1]. On the 
other hand, once this condition is fulfilled, B—L and CP 
violating interactions among ordinary and mirror parti- 
cles can generate baryon asymmetries in both sectors Q , 
naturally giving the relation il'^/fis ~ 5 between cos- 
mological fractions of the dark and visible matter '4] . 

The same B or L violating interactions that lead to pri- 
mordial baryogenesis can also induce mixing phenomena 
between the ordinary particles and their mirror partners. 
E.g. effective operator {1/M)lipl'(j)' (AL = 1) between the 
ordinary / mirror lepton and Higgs fields in the early uni- 
verse gives rise to an efficient lepto-baryogenesis mech- 
anism in both sectors, while at low energies it induces 
the mixing between ordinary (active) neutrinos Ve,p,,T and 
their mirror (sterile) partners v'^ ^ ^ [5^. 

Effective six-fermion interaction {1/M)^{udd){u'd'd') 
( AB = 1) involving the ordinary (u, d) and mirror (u' , d') 
quarks can provide an efficient baryogenesis mechanism 
at the scale M of few TeV which in addition can be 
testable at the LHC @. On the other hand, it induces 
the neutron— mirror neutron mixing, e{nn' -\- n'n). Since 
the masses of n and n' are exactly equal, in the vacuum 
conditions they must have a maximal mixing while the 
oscillation time r„„' = ^ (M/10 TeV)^ s in their rest 
frame can be much smaller than the neutron /3-decay time 
Td — 880 s; such a fast n — n' oscillation is not excluded 
experimentally 6] . This seems rather surprising as it re- 
gards to the baryon number violating process (AB = 1). 
The key moment is that for free neutrons n — n' oscilla- 
tion is suppressed by the Earth's magnetic field whereas 
it is ineffective for the neutrons bound in nuclei so that 
the nuclear destabilization limits become irrelevant. 



In the last years several experiments searched for n — n' 
oscillation via the magnetic field dependence of the neu- 
tron losses Note, however, that the lower bounds 
on r„„' reported by these experiments and adopted by 
the Particle Data Group [l3] (the strongest limit reads 
Tnn' > 414 s @) become invalid if the Earth possesses 
a reasonably large mirror magnetic field: the latter can- 
not be screened in the experiments and can drastically 
change the probability oi n — n' oscillation [ll| . On the 
other hand, if dark matter consists of mirror particles, 
it is also plausible that the solar system and the Earth 
itself may capture a significant amount of mirror matter 
due to feeble interactions between ordinary and mirror 
particles. The Earth's rotation gives rise to circular cur- 
rents in the captured mirror matter which may induce a 
mirror magnetic field up to several Gauss [ll| . 

In the presence of non-zero mirror field the proba- 
bility oi n — n' oscillation has a non-trivial dependence 
on the ordinary magnetic field *B and on its orienta- 
tion relative to 03' [ll|. In particular, for > 1 G, 
the strongest bound comes from measurements of the 
neutron losses in the Earth's magnetic field which yield 
w[s] x Q3'[G] > O.I [nl, [il]. Thus, for «B' - 10 G 
even t„„/ 10^^ s is allowed. Interestingly, such a fast 
n — n' oscillation is not in conflict with the primordial nu- 
cleosynthesis or the neutron star stability bounds [1, [HI ■ 
In addition, anomalous dependence of the neutron losses 
on the magnetic field orientation observed in the experi- 
ment [9] can be interpreted in terms oi n — n' oscillation 
with Tnn' of few seconds provided that 03' ~ O.I G [l^ . 

In this letter we show that fast n — n' oscillation can 
have intriguing implications for the propagation of ultra- 
high-energy cosmic rays (UHECR). It brings to signifi- 
cant modifications of the spectrum at i? > 10 EeV which 
can be proved in future with better accuracy. 

It is known that the cosmic microwave background 
(CMB) causes an abrupt end in the cosmic proton spec- 
trum, the so called Greisen-Zatsepin-Kuzmin (GZK) cut- 
off [iJl- The cutoff energy corresponds to the pion pho- 
toproduction threshold, Eqzk — nn^ml2e.y ~ 60 EeV, 
where m and m-^ are respectively the nucleon and pion 
masses and ~ ST ~ lO"'^ eV is an effective energy of 
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relic photons, T = 2.725 K being the CMB temperature. 
The mean free path (m.f.p.) of the proton, Ig ~ (o'(p7 
N-K)n-y)^^ cx T^^, strongly depends on its energy. One 
has Zs ~ 5 Mpc for E > 300 EeV but it sharply increases 
at lower energies becoming Ig > 100 Mpc at iJ < 60 EeV. 
In each pj Ntt scattering the super-GZK protons lose 
about 20% of their energy. 

The p7-scattering has two main pion production chan- 
nels, P7 — ^ p7r° and pj -> n7r+, with roughly comparable 
cross-sections. Conversion of the cosmic ray proton into 
the neutron does not influence the propagation length, 
since n^y Ntt scatterings, 71.7 — >■ mr^ {pTr~), have nearly 
the same cross sections as P7 — >■ Ntt ones. In addition the 
/3-decay n pei>e converts the neutron back to the pro- 
ton with practically the same energy. Up to ~ 0.5 ZeV 
the decay length 1^ = FcTd (F = E/m is Lorentz factor) 
is smaller than 717 — > Ntt scattering length. Hence, cos- 
mic ray carriers with E 3> Eqzk travel long distances 
transforming from protons to neutrons and back suffering 
significant energy losses which downgrade their energy to 
sub-GZK range. Yet, the baryon number in the cosmic 
ray propagation is conserved. 

In presence of n — n' oscillation the situation changes 
drastically: the produced neutron can now oscillate into 
a mirror one. If t„„' <C the oscillation length l^^r = 
cTTnn' is much smaller than and Is so that at these 
scales oscillations may be averaged. The n — n' transition 
probability reads [1, |Tl| 



^ ^ 2[l + (ra;w)^] 2 + g(S/100EeV)2' ^ ^ 

where lu = i|Ai„A«B| and A!B = <B - 03', Hn being the 
neutron magnetic magnetic moment and *8 and *8' being 
respectively the ordinary and mirror magnetic fields at 
the cosmological scales, or more precisely their transverse 
components. Factor q = 0.45 x (r„„'/ls)^ x (A<B/1 fG)^, 
shows the efficiency oi n — n' oscillation at £' ~ Egzk- 
Finally, /3-decay of mirror neutron n' p' e'v'^ converts 
a cosmic ray, being initially a proton, to a mirror proton. 

The latter can be converted to ordinary proton via in- 
verse chain of reactions: p'7' — nV scattering, n' — n 
transition and n peve- However, the propagation 
length of mirror protons is much larger than that of 
ordinary ones, 1'^ 3> h, as far as the temperature of 
mirror CMB is smaller than that of ordinary CMB, 
T'/T — X 1. Namely, the Big Bang nucleosynthe- 
sis imposes a robust upper bound x < 0.5 or so [2|, 
but the limits strengthen if one assumes that dark mat- 
ter consists entirely of mirror baryons. In this case the 
large scale structure and CMB power spectrum require 
X < 0.3, while yet stronger limits x < 0.2 {x < 0.1) arise 
by demanding that the Silk damping of mirror baryon 
perturbations does not prevent the formation of normal 
(dwarf) galaxies 

For the relic mirror photon number density and their 
average energy we have n!y = x^n^ and e'^ = xe-f. 
Thus m.f.p. of mirror cosmic rays drastically amplifies, 
I's ~ x~^ls, while the threshold energy of p'j' — > N'tt' 



becomes, E'q^^ ~ x~^Eqzk- So, the energy range 
Eqzk ^ E < i?Q2K ^'^ts for ordinary cosmic rays like 
a sink where they disappear - ordinary cosmic rays with 
E > Eqzk are converted to mirror ones, but the mirror 
ones may be converted (at much lower rate) to ordinary 
ones only at _B > -Eq2K- T^^^ dominant fraction of the 
cosmic rays produced in far distant sources must escape 
to the parallel sector via n — n' oscillation. However, if 
there are powerful mirror sources the ordinary UHECR 
flux may be increased at i? > £^q2k by the contribution 
from cosmic rays arriving from the mirror sector. 

In the presence of n — n' oscillation, evolution of 
the four UHECR number densities Ui = Ui{E,t), i — 
p,n,p' ,n' , in the expanding universe may be described 
by a system of coupled integro-differential equations 

°° ~ ~ (2) 

-R,{E,t)U,+T,j{E) jdEWjk{E,E,t)Uk{E,t), ^ ' 

E 

where Hit) is Hubble parameter. We assume that cos- 
mic rays sources are distributed homogeneously in space 
and their generation functions Qi(E, t) may have cosmo- 
logical evolution with time t [l3|. Here Wjk{E,E,t)) 
is the probability density for a nucleon k (N — p, n 
or N' = p',n') of energy E to transform via its pion- 
production scatterings off the respective CMB (7 or 7') 
into a nucleon j (again N or TV') with energy E. The 
matrix T.^j^E), with Tpp = Tp>p> = 1, r„„ = T„/„' = 
l-P{E), Tnn' = Tn'n = PiE) givcu by Eq. © and 
with other elements being zero, stands for transition 
probabilities due to n — n' oscillation: the neutron n 
produced in A'^7-scattering, N = p,n, promptly oscil- 
lates into n' with a probability P{E) and vice versa; 

Ri{E,t) = dE'£^Wji{E,E,t) stands for probabil- 
ity of a nucleon i with energy E to disappears from the 
energy range dE; factors /3p_p' (i?, t) take into account 
the p and p' energy losses due to e~^e~ pair produc- 
tion {f3n,n' = 0). The matrix Dij takes into account n 
and n' /3-decays transforming neutrons back into protons 
with practically the same energy; here Dpn = -Dp'n' = 1, 
Dnn = Dn'n' — and all other elements are zero. In 
the absence of n — n' oscillation the system ([2]) obviously 
splits into two independent sets of equations for ordinary 
and mirror cosmic rays. 

At all reasonable energies neutrons decay before their 
n7 Nil scattering off CMB, so that Id < k ^ I's- The 
relation 1^ <C l'^ holds very well for mirror neutrons, while 
for ordinary ones 1^ < Is is fulfilled only at < 0.5 ZeV. 
For this energy range the initial proton, after p7 — > n7r+ 
scattering, with following prompt n — n' oscillation and 
neutron decay, instantly transforms into a mirror proton, 
p ^ p' with probability PiE), and vice versa, p' — p, 
neglecting the propagation periods when a nucleon dwells 
in the mixed n — n' state. It should be noted that this 
approximation is perfect if the difference between pj- and 
n7-cross-sections is neglected. This allows to reduce the 
system of 4 equations ([5]) to a system of two equations 
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FIG. 1. Ratios of modified hy n ~ n' oscillation UHECR 
spectrum, J„„i{E), and "standard" spectrum, Jo{E) (no os- 
cillation), in the absence of mirror sources {Q'{E) = 0) and 
for different ratios of CMB temperatures x = T'/T. Source 
generation spectrum is taken as Q{E) oc -E~^'0(iJmax — E). 



describing evolution of just p and p' . 

To facilitate the calculations, the latter system of two 
integro-difFerential equations was reduced to an analog 
of a set of two coupled Fokker-Plank type differential 
equations by expanding the kernels of integrals in ^ in 
series at i? = i? up to second derivatives. This method 
was proved to be valid by comparison with Monte Carlo 
simulations in the case of propagation of ordinary pro- 
tons [l^. Moreover, we assume g ^ 1 in Eq. ([1]), so that 
P{E) = 1/2. This holds for A«B < 1 fG. The data con- 
cerning the cosmological magnetic fields are controversial 
and there are some hints that they may be larger. How- 
ever, at large scales ordinary and mirror magnetic fields 
can be strongly correlated, so that their difference is van- 
ishing 6]. There is also a possibility of the resonance 
MSW like transitions if magnetic fields have turbulent 
structure at scales less than 1 Mpc. 

The results of our calculations for different sets of pa- 
rameters are shown in Fig. [T] and Fig. [5] Let us first 
discuss the case when there are no mirror cosmic rays 
sources, Q' = 0. Then cosmic rays with E > Eqzk pro- 
duced in distant extragalactic sources not only lose their 
energy during propagation, but also degrade in number 
owing to n — n' transition. This makes the cutoff in 
the cosmic ray spectrum to begin at energies lower than 
Eqzk — 60 EeV. Now the cutoff relates also to the non- 
conservation of the baryon number B: the most part of 
initial protons transforms into mirror protons, thus get- 
ting invisible. On the other hand, an imaginary mirror 
observer would detect a flux of mirror protons with ener- 
gies E > Eqzk originated in our world. In other words, 
the baryon numbers B and B' are not conserved individ- 
ually, but the sum B + B' must be conserved. 

In the Fig. [1] we show the modification factor ri{E) ~ 
Jnn'{E)/ Jq{E), defined as the ratio of UHECR spectra 
calculated with and without n — n' oscillation, in the 



FIG. 2. The same as on Fig. [T]but in the presence of mir- 
ror sources, Q'{E) with 7g = 2.5, different ratios of source 
intensities, Q'{E)/Q{E) = 0,1,2,5, and with T'/T < 0.3. 



absence of mirror sources, Q' = 0. Actually, it has an 
almost universal shape just weakly depending both on 
the source generation function index 7^ and on the ratio 
of CMB temperatures x = T'/T in two sectors. Re- 
ally, the spectrum cutoff with oscillation starts earlier, 
at i? ^ 10 EeV, while the maximal difference from the 
GZK prediction (about 40%) is reached around 100 EeV. 

Fig.Oshows the same spectral modification factor ri[E) 
in the presence of the mirror sources, for different values 
of Q'{E)/Q{E). In fact, if mirror baryons constitute dark 
matter, i.e. il'^/fls, one can expect that there are also 
mirror cosmic rays. In fact, in quasars and AGN the cen- 
tral black holes can be a potential sites for acceleration of 
both ordinary and mirror protons. Therefore, it is natu- 
ral to assume that the ratio Q'{E)/Q{E) is the same as 
the ratio of the relative matter fractions, fl'^/fls- Let 
us recall, that n'g/flB ^ 5 can be naturally achieved 
in the joint ordinary-mirror baryogenesis mechanism Q 
provided that x < 0.3 or so 

In this case, in spite of larger m.f.p., the mirror cos- 
mic rays can partially move to the ordinary sector and 
substantially increase the UHECR flux at energies above 
E = 100 EeV. The value of the turning point depends 
on the parameter x = T'/T as well as on the strength 
of mirror sources, Q'{E)/Q{E). Note however that the 
existence and position of the pre-GZK cutoff is robust: 
the shape of the r]{E) at i? = 10 — 100 EeV practically 
does not depend on the strength of the mirror sources. 

Unfortunately, the data statistics achieved at the op- 
erating installations is not yet enough to determine the 
exact shape of the spectrum above 100 EeV. However, 
the earlier end of the cosmic ray spectrum seems to be 
indicated by the data of PAO ^Ifi] . A controversy between 
the Auger and HiRes [i3 - TA [H] data concerning the 
shape and chemical composition of the UHECR spectrum 
is not yet resolved (see e.g. [l^). Auger sees the cutoff 
starting from ~ 25 EeV, a factor of 2 lower than Eqzk 
and claims to see here heavy nuclei. The estimated sys- 
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tematic uncertainty of 22% is not sufficient for adjusting 
the cutoff position. The HiRes data indicate cosmic rays 
to be protons with cutoff at Eqzk, but their statistic is 
lower. It is worth to note, however, that the issue of 
the composition determining at such high energies is a 
subject of theoretical models rather than of experimen- 
tally measured cross-section and multiplicities (see e.g. 
[20|). The proposed n — n' oscillation model mitigates 
the controversy in the data of two experiments putting 
the beginning of the cutoff in the middle, though extra- 
galactic protons are assumed to be cosmic ray carriers. 
It should be noted, that in case the UHECR are mostly 
heavy nuclei, the neutrons produced in their photodisin- 
tegration on CMB will be also lost as they escape to the 
parallel world. 

Prompt neutron oscillations may provide correlation of 
cosmic rays with E > 100 EeV with distant sources (e.g. 
BL Lacs [2l[) as an indication of the UHECR transport 
in the parallel world. It is interesting to remark, that 
cosmic ray fraction which comes from the parallel sector, 
can be composed by anti-protons rather than by protons. 

Another immediate consequence of the baryon losses 
due n — n' transitions is a strong suppression of the 
cosmogenic Berezinsky-Zatsepin (BZ) [2^ neutrino flux 
mostly produced via P7 — > mr^ scattering with following 
7r+ — >■ /x+i'^ and /i+ e'^v^v^ decays. This conclusion 
remains pessimistic even if ordinary and mirror neutrinos 
also have non-zero mixing Q. Really, since T' <T, mir- 
ror cosmic rays have very large m.f.p., l'^ ~ x~^ls (e.g. 



I'g ~ 600 Mpc for X = 0.2 versus Ig ^ 5 Mpc), and thus 
suffer much less scatterings than ordinary ones. There- 
fore the diffuse cosmogenic neutrino flux may turn out 
to be much lower than expected [23]. Even the giant 
ICECUBE 24] with its control over 1 km'^ of ice may be 
insufficient to detect this flux. 

The electromagnetic cascades originated from — > 
p7r° channel with subsequent tt" — > 27 decay and 
77CMB ^ e+e~ will be also suppressed, since most 
amount of the ordinary super-GZK cosmic rays escape 
to the parallel world just after few proton scatterings 
off CMB. This reduces the restrictions imposed on the 
UHECR models by the diffuse extragalactic gamma-ray 
flux measured by Fermi-LAT ^ at E > 100 GeV 

To conclude, the interesting effect of the fast baryon 
number violation via n — n' oscillation may be detected in 
the laboratory neutron disappearence and regeneration 
experiments, while the underlying TeV-scale physics can 
be tested at the LHC ^] . Here we show that this attrac- 
tive concept may complementary bring to the testable 
modiflcations of the cosmic ray spectrum. Due to the 
baryon non-conservation during the UHECR propaga- 
tion the cutoff of the spectrum shifts to lower energies 
and becomes signiflcantly steeper. On the other hand, 
in the presence of powerful mirror sources the observed 
spectrum at highest energies may get even higher as com- 
pared to the 'no oscillation' case. 
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